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Abstract. In 2008, the blazar Mrk421 entered in a very active phase and was one of the 
brightest sources in the sky at TeV energies, showing strong and frequent flaring. We 
searched for y-ray emission at energies E>0.8 TeV during the whole 2008 with the ARGO- 
YBJ experiment, a full coverage air shower detector located at Yangbajing (4300 m a.s.l., 
Tibet, RR. China). The observed signal is not constant and in correlation with X-ray mea- 
surements. The average emission, during the active period of the source, was about twice 
the Crab Nebula level, with an integral flux of (4.9±2.0)xl0~ n y crrT 2 s _1 for E r >l TeV. 
This paper concentrates on 2008 June when the Mrk421 flaring activity has been studied 
from optical to 100 MeV gamma rays, and only partially up to TeV energies, since the 
moonlight hampered the Cherenkov telescope observations after 8 June. Our data complete 
these observations, with the detection of a second flare of intensity of about 7 Crab units 
on June 11-13, with a statistical significance of 4.2 standard deviations. The observed flux 
is consistent with a prediction made in the framework of the Synchrotron Self-Compton 
model, in which the flare is caused by a rapid acceleration of leptons in the jet. 
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1. Introduction 

Mrk421 (z=0.031) is the brightest BL Lac 
object and the first extra galactic source dis - 
covered at TeV energies dPunch et al. 1992I) . 
where dramatic variability has been observed 
with flux increasing by more than a fac - 
tor 50 in about one hour dGaidos et al. 1 996). 
Like most blazars, its spectral energy distri- 
bution (SED) shows two smooth broadband 
components, the first one peaking in the soft 
to medium X-ray range and the second one 
extending to the GeV/TeV energies. X-rays 
are generally attributed to synchrotron radi- 
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ation from high energy electrons, while the 
origin of the y-ray emission is more un- 
certain. Possibilities include inverse Compton 
scattering of synchrotron (Synchrotron Self- 
Compton, SSC) or ambient photons (external 
Compton, EC) off a single electron population 
(Samb runaet al. 19961: iGhisellini etal. 19981 
iFossati et al. 20081) . Nevertheless, an alterna- 
tive hadronic m odel (y-rays from proton syn- 
chrotron dMucke et al. 20031) ) for Mrk421 is 
not ruled out yet. 

Multiwavelength observations are the 
key for understanding the blazar phe- 
nomenon. Indeed, it is generally true that 
hadronic models are in trouble to re- 
produce the observed highly correlated 
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X-ray/TeV variability, t hat strongly sup- 
ports the SSC models dFossati et al. 2008t 
iFidelis & Iakubovskvi 20081: IWagner 20081) . 

The ARGO-YBJ experiment is presently 
the only wide-field of view y-ray telescope 
able to detect AGN TeV flaring activity on 
a few day period. In this paper we report 
on the monitoring of Mrk42 1 performed with 
ARGO-YBJ in 2008. We will give special at- 
tention to the 2008 June flaring activity be- 
cause an extraordinary set of simultaneous 
measurements covering 12 decades of energy, 
from optical to TeV gamm a rays, is available 
dDon narumma et al. 2009). 



2. The ARGO-YBJ experiment 

The ARGO-YBJ detector, located at the 
YangBaJing Cosmic Ray Laboratory (4300 m 
a.s.l., Tibet, P.R. China), is the only exper- 
iment exploiting the full coverage approach 
at very high altitude. The detector is consti- 
tuted by a central carpet ~74x78 m 2 , made 
of a single layer of Resistive Plate Chambers 
(RPCs) with ~92% of active area, enclosed 
by a partially instrumented guard ring that ex- 
tends the detector surface up to ~100xll0 
m 2 . The apparatus has a modular structure, 
the basic data acquisition element being a 
cluster (5.72x7.64 m 2 ), divided into 12 RPCs 
(2.80x1.25 m 2 each). Each chamber is read 
by 80 strips of 7x62 cm 2 (the space pixel), 
logically organized in 10 independent pads of 
56x62 cm 2 representing the time pixel of the 
detector. The RPCs are operated in streamer 
mode with a standard gas mixture (Argon 15%, 
Isobutane 10%, TetraFluoroEthane 75%), the 
High Voltage settled at 7.2 kV ensures an over - 
all efficiency of about 96% (lAielli et al. 20061) . 
The full detector is composed of 153 clus- 
ters for a total active surface of ~6700 m 2 . 
All events giving a number of fired pads 
Nparf > N, r ,-g in the central carpet within a 
time window of 420 ns are recorded. The spa- 
tial coordinates and the time of any fired pad 
are used to reconstruct the position of the 
shower core and the ar rival direction of the pri- 
mary, as described in (|Pi Sciascio et al. 20071 
iDi Sciascio et al. 20081) . 



Since 2007 November the full detector is 
in stable data taking at the multiplicity trigger 
threshold N rn? >20 with a duty cycle ~ 90%: 
the trigger rate is about 3.6 kHz. 

3. Data analysis 

The dataset for the analysis of Mrk421 pre- 
sented in this paper contains all showers with 
N pad >40 and zenith angle 9 <40°. 

A 20° x 20° sky map in celestial coordi- 
nates (right ascension and declination) with 
0.1°x0.1° bin size, centered on the source 
location, is filled with the detected events. 
The background is evaluated with the t ime 
swapping method (lAlexandreas et al. 19921) . N 
"fake " events are generated for each detected 
one, by replacing the measured arrival time 
with new ones. These times are randomly se- 
lected within a 3 hours wide buffer of recorded 
data. Swapping the time means swapping the 
right ascension, keeping unchanged the decli- 
nation. A new sky map (background map) is 
built by using 10 such fake events for each real 
one, so that the statistical error on the back- 
ground can be kept small enough. 

To maximize the signal to noise ratio, the 
bins are then grouped over a circular area of 
radius iff, i.e. every bin is filled with the con- 
tent of all the surrounding bins whose center 
is closer than if/ from its own center. When 
the Point Spread Function of the detector is 
a Gaussian with r.m.s. <x, the opening an- 
gle if/ containing 71.5% of the events maxi- 
mizes the signal to background ratio for a point 
source with a uniform background, and it is 
equal to 1.58-cr. The values of iff are 1.9°, 0.9° 
and 0.5° respectively for N pad >40, 100 and 
300, in agreement with M onteCarlo simula- 
tions dVernetto et al. 2009al) . 

Finally, the integrated background map is 
subtracted from the corresponding integrated 
event map, thus obtaining the "source map". 
For each bin of this latter map, the significance 
with respect to the background is calculated. 

This analysis procedure has been tested 
with the Crab Nebula, the standard candle for 
VHE astronomy. At the Yangbajing latitude the 
Crab culminates at zenith angle 9 cu im - 8.1° 
and it is observable every day for 5.8 hours 
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Fig. 1. The Crab Nebula energy spectrum measured 
in 2008 by ARGO-YBJ compared with the results of 
some other detectors. 
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Fig. 2. The Mrk421 energy spectrum measured by 
ARGO-YBJ from 2008 day 41 to 180, when the 
source was in active state. The line represents the 
best fit to the data. 



with a zenith angle 6 <40°. The Crab Nebula 
has been observed from 2007 November to 
2009 March, for a total of 424 days on-source, 
obtaining a signal with a statistical significance 
of 7 standard deviations for N pai j >40. The cor- 
responding photon median energy is 1 . 1 TeV. 
In this analysis no events selection or y/hadron 
discrimination algorithm has been applied. The 
average number of gamma rays detected per 
day in the observational window centered on 
the source position is 155+25. 

By assuming a differential spectrum dN/dE 
= K E a , we simulated a source in the sky fol- 
lowing the diurnal path of the Crab, and eval- 
uated the number of events expected in the 3 
N pad bins 40-99, 100-299 and >300, for differ- 
ent values of K and a (10 -11 <K<10 10 photons 
cirT 2 s" 1 TeV- 1 and 1.5< a <3.5). The multi- 
plicity bins correspond respectively to median 
energies 0.8, 1.8 and 5.2 TeV. Comparing the 
expected values with the observed ones we ob- 
tained the spectrum that best fits to data 

dN/dE = (3.7 + 0.8) x io- n £- 2 - 67±0 - 25 (l) 

in fair agreement with other experiments (see 
Fig. [Q- 



4. Mrk421 analysis 

The same analysis has been performed for 
Mrk42 1 . This source culminates at the ARGO- 
YBJ location at zenith angle 6 cu i,„ = 8.1°, and 
it is observable every day for 6.38 hours with a 
zenith angle 6 < 40°. 

Here the Mrk421 data collected from 2007 
day 3 1 1 to 2008 day 366 are presented. Using 
the same method adopted for the Crab Nebula, 
we evaluated the Mrk42 1 spectrum from 2008 
day 41 to day 180, where the X-ray flux 
showed the most intense flares. In this period 
(755 observation hours) the observed signal 
had a statistical significance of 6.1 standard 
deviations. The differential spectrum (photons 
cirT 2 s -1 TeV -1 ) that best fits to data is (see 
Fig. |2) 

dN/dE = (7.46+1.70)xl0 _11 £ _2 - 51±0 - 29 e _T(£) (2) 

where the exponential factor e~ T(E> takes into 
account the absorption of gamma rays in the 
Extragalac tic Background Li ght, with r(E) 
given in (iPrimack et al. 20051) . The integral 
flux above 1 TeV is (4.9+2.0)xl0" n photons 
cirT 2 s _1 , about twice the Crab Nebula flux. 
The median energies corresponding to the 3 
N pad bins 40-99, 100-299 and >300 are re- 
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Fig. 3. Plot of the RXTE/ASM X-ray and ARGO- 
YBJ TeV gamma-ray count rates of Mrk421. The 
line shows the best linear fit. 



spectively: E = 0.84™^ , 1.74™-*| and 4.2™ ™ 
TeV. 

The observed gamma ray rate appears to 
be correlated with the X-ray rate measured 
by the All Sky Monitor detector aboard the 
RXTE satellite in the 1.5-12 keV energy range 
dRXTE web siteh . as can be seen in Fig. [3] 
where the X-ray and ARGO-YBJ TeV gamma- 
ray count rates are shown for all the simultane- 
ous measurements in 2008. The excess events 
observed by ARGO-YBJ, averaged over 10 
days, refer to showers with N pac i > 100. A posi- 
tive correlation between the rates seems appar- 
ent: the correlation coefficient is 0.64. 

Since Mrk421 is known to vary on different 
time scales, the source has been studied dur - 
ing 1, 10 and 30 days dVernetto et al. 2009bl) . 
For this analysis we have considered the data 
taken in the period 2007 day 311 - 2009 day 
89. No excess has been observed on a daily 
scale. Concerning the 10 scale, we observed an 
excess at 4.6 s.d. in the time interval 2008 days 
161 - 170, during a strong X-ray flare. Looking 
for 30 days excesses, the search has been car- 
ried out by shifting the 30 days intervals in 
steps of 10 days. We found several excesses 
from Mrk421 with significances between 4 and 
5 s.d., in particular in the intervals: 2008 days 



1 - 30, 71 - 100, 81 - 110, 91 - 120, 141 - 170, 
when several X-ray flares have been observed. 

4.1. The June 2008 flares 

A set of simultaneous measurements cover- 
ing 12 decades of energy, from optical to 
TeV gamma rays, was performed during the 
strong flaring activity in the first half of June 
2008 by different detectors: WEBT (optical 
R-band), SWIFT (UV, soft and hard X-rays), 
RXTE/ASM (soft X-rays), AGILE (hard X- 
rays and gamma rays) and the Cherenkov tele- 
scopes VERITAS and MAGIC (VHE gamma 
rays) (Donnarumm a et al. 2009). 

In this period two flaring episodes were 
reported, the first one on June 3-8, ob- 
served from optical to TeV gamma rays, the 
second one, larger and harder, on June 9 
- 15, observed from optical to 100 MeV 
ga mma rays. Using this mult i-frequency data, 
in dDonnarumma et al. 20091) the authors de- 
rived the SED for June 6, that shows the typical 
two humps shape, in agreement with the SSC 
model. According to the authors, the second 
hump intensity (that reached a flux of about 3.5 
Crab units at energy E >400 GeV) seems to in- 
dicate that the variability is due to the harden- 
ing/softening of the electron spectrum, and not 
to the increase/decrease of the electron den- 
sity. Their model predicts for the second flare a 
VHE flux about a factor 2 larger with respect 
to the first one. Unfortunately there were no 
VHE data included in their multi-wavelength 
analysis after June 8 because the moonlight 
hampered the Cherenkov telescopes measure- 
ments. 

The VHE observation was actually made 
by the ARGO-YBJ experiment, that since 2007 
November is performing a continuous mon- 
itoring of Mrk421. Fig. H] shows the rate of 
events with N pad >100 observed by ARGO- 
YBJ during the first half of June, averaged over 
3 days, compared with the X-ray flux measured 
by RXTE/ASM. 

A second flare has been detected with a 
statistical significance of 3.2 standard devia- 
tions during the interval 11-13 June. The sig- 
nificance increas es to 4.2 s.d. with a suitable 
events selection dVernetto et al. 2009ah . Fig. [5] 
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Fig. 4. Upper panel: rate of events with N piK j >100 
observed by ARGO-YBJ as a function of time from 
May 31 00:00 UT to June 19 00:00 UT. Each bin 
contains the rate averaged over the 3 days interval 
centered on that bin. Lower panel: daily counting 
rate of RXTE/ASM. 



shows the 6°x6° sky map around the source 
position in these 3 days, after applying the 
event selection. For every 0.1°x0.1° bin, the 
map gives the value of the statistical signifi- 
cance of the excess of events inside the circu- 
lar window of radius 0.9° centered on that bin. 
Donnarumma et al. evaluate a theoretical SED 
curve for the first flare fitting the observations 
made on June 6 from optical up to VHE en- 
ergies. During the second flare, they reported 
a higher photon flux from soft X-rays to 100 
MeV gamma rays. From these data they pre- 
dict a flux at E>1 TeV of 1.45-10" 10 photons 
cirT 2 s _1 corresponding to about 7 Crab units, 
and they model a SED curve with the Inverse 
Compton hump slightly shifted towards higher 
energies. 

Fig. |5] shows the event rate observed by 
ARGO-YBJ (in 18.2 hours of measurement) 
compared with the rate expected from a source 
spectrum given by the theoretical SED, for any 
Npad interval. The agreement is good. A simi- 
lar analysis is made for the first flare, integrat- 
ing our data on June 5-7 (17.9 hours of mea- 
surement). The observed signal has a signifi- 
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Fig. 5. Event rates observed by ARGO-YBJ as a 
function of the minimum pad multiplicity on June 
5-7 and June 11-13 (respectively triangles and cir- 
cles) compared with expected rates according to the 
Donnarumma et al. model for the same two periods 
(respectively dashed and solid lines). The inset rep- 
resents the sky map around the Mrk421 position on 
June 11-13, obtained for events with N p(lrf >100. The 
color scale shows the significance of the signal in 
standard deviations. The circle represents the obser- 
vational window of radius 0.9°. 



cance of ~2 standard deviations, increased to 
3 using the data selection. The event rate ob- 
tained as a function of the minimum pad mul- 
tiplicity, is consistent with the one predicted by 
the SED of Donnarumma et al. for June 6. 

Finally we estimate the energy spectrum 
for the second flare (the marginal signifi- 
cance of the first flare doesn't allow the spec- 
trum evaluation). Assuming a source spec- 
trum given by the theoretical SED, the me- 
dian energy of the events in the 3 pad mul- 
tiplicity bins (40 - 99, 100 - 299 and >300) 
are respectively 0.9, 1.4 and 2.4 TeV. Fig. 
[6] gives the measured SED for the 3 energy 
points ((3.38+2.03)-10- 10 , (3.55+1. 37)-10- 10 
and (2.49+1. 63)- 10~ ll) erg _1 cm 2 s), together 
with all the measurements in the optical-TeV 
range, and the theoretical SED for the two 
flares. The measurements appear in fair agree- 
ment with the expected emission. 
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Fig. 6. SED measured by ARGO-YBJ on June 
11-13 (large blue filled circles) together with the 
data of other experiments, obtained during the first 
flare (open circles), and the second one (filled cir- 
cles). The curves represe nt the SEDs modeled in 
dDonnarumma et al. 2009h for the first flare (dashed 
line) and the second one (solid line). The inset 
shows a zoom on the ARGO-YBJ data. 



5. Conclusions 

In summary, Mrk421 has been continuously 
monitored with ARGO-YBJ during 2008, 
showing a VHE flux twice the Crab Nebula 
level from day 41 to 180, when the source was 
in active phase, and decreasing afterwards. 

ARGO-YBJ observed a flare on 5-7 June, 
with a flux about 3.5 Crab units, in agreement 
with the VERITAS observation. A second flare 
has been discovered with a statistical signif- 
icance of 4.2 s.d., during the interval 11-13 
June, with a flux about 7 Crab units. 

ARGO-YBJ measured the spectra of 
Mrk421 above 0.8 TeV during the second flare 
completing a multiwavelength campaign from 
optical to TeV energies. 

For the first time an air shower array was 
able to detect gamma-ray flaring activity at 
sub-TeV energies on a few days period. 
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